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Abstract 

A two-sided exit problem is solved for a difference of a compound Pois- 
son process and a compound renewal process. More precisely, the Laplace 
transforms of the joint distribution of the first exit time, the value of the 
overshoot and the value of a linear component at this instant are found. 
Further, we study the process reflected in its supremum. We determine the 
' main two-boundary characteristics of the process reflected in its supremum. 

\ These results are then applied for studying the M'^jC'^llIB system. We 

■ derive the distribution of a busy period and the numbers of customers in 
\ the system in transient and stationary regimes. The advantage is that these 

■ results are in a closed form, in terms of resolvent sequences of the process. 



Introduction Queueing systems with batch arrivals and finite buffer 
have wide applications in the performance evaluation, telecommunications and 
manufacturing systems. One of the crucial performance issues of the single-server 



^ . queue with finite buffer room is losses, namely, customers (packets, cells, jobs) 



that are not allowed to enter the system due to the buffer overflow. This issue is 
especially important in the analysis of telecommunication networks. Motivated 
by this fact, we derived the most important performance measures of queueing 
systems of this type. More precisely, we considered the M''\G^\l\B queuemg 
system with finite buffer and its modification. We consider partial rejection, 
meaning that if an overfiow of buffer occurs due to the arrival of a batch of 
customers, the amount of work brought by this batch is only partially admitted 
to the buffer, up to the limit of the free buffer space just before the arrival. The 
rest is rejected and therefore, is lost. 

Evolution of the number of customers in such systems is described by a process 
with two refiecting boundaries. In general case this process is a difference of two 
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compound renewal processes. Reflections from the upper boundary are generated 
by the supremum (infimum) of the process. Reflections from the lower boundary 
govern the server's behavior. In general, such processes are not Markovians, 
but by adding a complementary linear component (in some literature called age 
process), we obtain a Markov process, which describes functioning of the queueing 
system. Studying main characteristics of the system results to the investigating 
the two-boundary functionals of the governing process. We applied the solutions 
of the two-sided exit problem for the governing process (see [25j and [29\ for 
the methodology) to obtain performance measures of interest. For the queueing 
system of M^|G^|l|i?, G^|M^|1|B (see [28]) type the governing process is the 
difference of the compound Poisson process and the compound renewal process 
complemented with the age process. The main result of this paper is the closed 
form formulae for the Laplace transforms of the busy period, time of the first loss 
and the number of customers in the system at arbitrary time. 

First passage times of the level by Levy processes in context of queues were 
considered in [12j, where the explicit characterization of the Laplace transform 
of the busy period distribution was found for a finite capacity M|G|1 queue, see 
also [H], and [6]. In regard to finding the buffer overflow time we refer to 
[2], where arrivals are modeled by a Markov modulated Poisson process (MMPP) 
and service time is exponential, and also to pjj, where BAIAP\G\l\b queue was 
considered. Previous works on the overflow period were concentrated on simple 
Poisson arrivals [7], [lOj . or renewal arrivals |14) |42j . Picard and Lefevre |45] 
found the distribution of the first crossing time of a Poisson process and renewal 
process in terms of polynomials of Abel-Gontcharoff types. 

In recent years there has been a great interest in analyzing various queueing 
models with MAP (Markov Arrival Process) as input process or MSP (Markov 
service process). MAP is used to represent correlated traffic arising in modern 
telecommunication networks. In systems with Markov arrival or service processes 
(MAP, BMAP, or BMSP) and their modifications, it is common to use the sup- 
plementary variable methods and/or embedded Markov chains. For the method 
of supplementary variable we refer, for instance, to [15j, where GI\MSP\1 queue 
with finite as well as infinite buffer was analyzed. Embedded Markov chains tech- 
niques were used by [IT]. De Boer et.al. [7j studied stationary distribution of 
the remaining service time upon reaching some target level in an M|G|1 sys- 
tem. The asymptotic analysis of G\MSP\l\r queue has been carried out by [8]. 
Banik et.al. [3] found steady state distribution for a finite-buffer single-server 
queue GI\BAISP\1\N with renewal input. Random size batch service queueing 
models were subject of study in [13], [9] (stationary analysis of GI\BMSP\1 
queue) and [30] (asymptotic behavior of the loss probability). 

Hence, the majority of recent literature is devoted mainly to the queue size 
and workload, most of the times in the steady state case. However, recently it was 
shown that steady-state parameters do not reflect the reality. A detailed discus- 
sion of the drawbacks of steady-state parameters in telecommunication networks 
may be found in [48]. This remark emphasizes importance of studying main per- 
formance measures in transient regime, which is the topic of this article. The 
main two-boundary characteristic of the governing process is the joint distribu- 
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tion of {x,L,T} , i.e. of the first exit time from the interval, the value of the 
overshoot and the value of the linear component at this instant. For the overview 
of the existing results on the two-boundary problems we refer to [25] . And here 
we only mention several authors who contributed a lot in the development of this 
area. Starting from Kemperman (1963), Takacs (1966), Emery [TU], Pecherskii 
[33] . Suprun, Shurenkov ([50], [51]), Lambert [35], Doney [16], Avram, Kypri- 
anou, Pistorious [3], Pistorious [46], Kyprianou, Palmowsky [34], and Kadankov, 
Kadankova ( [23 , [27] ) studied one- and two-boundary characteristics for different 
classes of stochastic processes. 

The Laplace transforms of the joint distribution of the first exit time and 
the value of the overshoot at this time instant for general Levy processes and 
random walks have been determined in [29]. The Laplace transforms of this 
joint distribution were found in terms of the Laplace transforms of the one- 
boundary characteristics of the process. This method for Levy processes and 
random walks [29| was then applied for other classes of stochastic processes, such 
as the difference of compound renewal processes [27j, and semi-Markov random 
walks with linear drift |24j . 

The rest of the article is structured as follows. In Section 1 we introduce the 
process, necessary notation and consider the one-boundary characteristics of the 
process. The two-sided problem is solved in Section 2. In this section we also 
prove the weak convergence of the joint distribution of the supremum, infimum 
and the value of the process to the corresponding distribution of the symmetric 
Wiener process. Section 3 deals with the reflected processes. We consider the first 
passage of the lower boundary, distribution of the increments of the process and 
its asymptotic behavior. Finally, in Section 4 we apply the results obtained in 
the previous sections for studying several characteristics of the queueing system 
M^|G^|1|B, such as busy period, time of the first loss of a customer and the 
number of customers in the system in transient and stationary regimes. 

1 Preliminaries and definitions 

Let X, 5 G N = {1,2,...} be positive independent integer random variables 
and let rj £ (0, oo) be a positive random variable independent of x, 6 with 
the distribution function F{x) = P [// < x] , x > 0. We will assume that Ex, 
EJ, E?7 < oo. Introduce the sequences {r],r]'^}, {x, >4}, {6,5'^}, n G N of 
independent identically distributed (inside of each sequence) variables and define 
the monotone sequences 

r]o{x) = 0, r]i{x) = 7]^, r]n+i{x) = Vx + Vi ^ ??n, n G N, (1) 

xo = 0, x„ = x^ H h x^; (5o = 0, 5n = 6[ -\ h 5^; n£ N, 

where r/^. G (0, oo) is a random variable with the following distribution function 
F^{u) = P [r]^ <u\ = [F{x + u)- F{x)]{l - F{x))~'^ n > 0. 

Denote by {T^{t)}t>o £ '^'^ = {0, 1, . . . } a compound Poisson process with the 
generating function of the form 

^Q^{t) ^ ^tm^ k(e) = niEO'' -1) , \9\<i, 
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where fi > is the intensity of the jumps and x is a jump size. For all t > 
define a renewal process generated by the random sequence {r]nix)}n&z+ ^ 
follows 

N^{t) = max{n G Z+ : r?„(x) < t} G Z+, x> 0. 
Introduce a right-continuous step process for all x > 

Z),(t) =7r(t)-(5^^(t) GZ = {0,±1,...}, t>0; D,{0) = 0. (2) 

Note, that inter-arrival times of the positive jumps are exponentially distributed 
with parameter n, the positive jumps themselves are of a random size x, and 
there occur negative jumps of size 5'^ at time instants r]n{x), n G N. We will 
call the process {Dx{t)}t>o a difference of the compound Poisson process and a 
compound renewal process. Observe, that this process is not a Markov process 
in general. For all t > introduce a right-continuous linear component 

r? + (t) = | J + GM+ = [0,oo), x>0. (3) 

The process {il:^i't)}t>o increases linearly on the intervals [r]n{x), r]n+i{x)), nG 
Z+, it is killed to zero at the points r]n{x), n G N, and the value of the pro- 
cess at the instant to ^ Vx is equal to the time elapsed from the moment of 
the last negative jump of the process ([2]) till Iq. We will call the process ([3]) a 
linear component (sometimes referred to as the age process). By adding this lin- 
ear component to the process {Dx{t)}^yQ we obtain a right-continuous Markov 
process 

{Xt}t>o = {Dx{t),r] + it)}^^^€ZxR+, Xo = {0,x}, x > 0, (4) 

which governs the process {Dx{t)}t>o- The process defined in (jH) is a Markov 
process. Note, that it is homogeneous with respect to the first component [18j. 
This means that if Xt^ = {k, u}, k & 1^, n > 0, then the evolution of the process 
{Xt}t>to in the sequel does not depend on the value k of the first component, 
and the first positive jump of the process {Dx{t)}t>to (which is distributed as 
x) will occur after an exponential period of time with parameter fi. The first 
negative jump of the process {Dx{t)}t>to (which is distributed as 5) will take 
place after elapsing of time This fact will be used constantly when setting 
up the equations. 

Here and in the sequel we assume that the random variable 5 G N is geomet- 
rically distributed with parameter AG [0, 1) : 

P[J = n] = (1 - A)A"-\ nGN, BO^ = 9^—^, \e\<l. 

1 — xo 

This assumption means that the process {Dx{t)}t>o has geometrically distributed 
negative jumps at time instants {??n(a^)}neN- Throughout the article we will use 
the following notation 6 ~ ge{X). In this case it is possible to obtain closed form 
solutions for the one and the two-sided boundary problems. Our task now is 
to determine the Laplace transforms of the joint distributions of the upper and 
lower one-boundary functionals of the process {Xt}t>o- In the sequel we will use 
the following result. 
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Lemma 1. Let f{s) = Ee Then for s > the equation 

0-X = (l-X)f(s-k{9)), \9\<l (5) 

has a unique solution c(s) inside the circle \6\ < 1. This solution is positive 
and c{s) G (A, 1). If B[>i],B[r]] < oo, p = - A)E[>ir] E[r/], then for p > 1, 
lim c(s) = c G (A, 1); and for p < 1, lim c(s) = 1. 

A detailed proof of an analogous proposition for semi-continuous random 
walks can be found in the monograph of Spitzer [17] . The reasoning in that proof 
can be applied to the equation (0) as well (see also Lemma 1 [27J). 

Let Xo = {0,x}, xGM+, k€Z+. Define 

rfc(x) = inf{t : D^{t) < -k}, Tk{x) = -D^{Tk{x)) - k, inf{0} = oo, 

i.e. the first overshoot time of the negative level —k by the process {Dx{t)}t>o- 
We will use the convention that on the event {Tfc(x) = oo} Tk{x) = oo. Denote 
5Sfc(x) = {rfe(x) < oo}. 



A(x,m,s) =E e-^"'=(");Tfc(x) = m,«fc(x 



m G N. 



The Laplace transforms of the joint distribution of the lower one-boundary func- 
tionals are determined by means of the following lemma. 

Lemma 2 ([25j). Let {-Dx(t)}t>o the difference of the compound Poisson 
process and the compound renewal process, 5 ~ ge{\). Then 

(i) the Laplace transform of the joint distribution of {Tfc(x), Tyfc(x)}, k G Z"^, 
X > satisfies the following equality for s > 0, m G N 

fk{x,m,s) = Us - k{c{s)))c{s)^ (1 - (6) 

where c{s) G (A, 1) is the unique solution of the equation ^ inside the 
circle \e\ < 1, Us) = Ee-'^''^ f{s) = Ee"^" = Us); 

(ii) if p > 1, then P[Tk{x) < oo] = fx{—k{c))c^ < 1, and Tk{x) for all 
k G X > is a defective random variable; if p < 1, then P[Tfc(x) < 
oo] = 1, and Tk{x) is a proper variable for all k G Z"^, x > . 

Observe that for all k £ Z the value of the overshoot Tf^ix) does not depend 
on Tfc(x) and it is geometrically distributed Tfc(x) ~ ge{\). 

We now introduce a sequence which will be used to obtain the results in the 
sequel. The idea to employ this sequence for semi-continuous random walks and 
semi-continuous Levy processes is due to Takacs [53]. Since the function 



Us - k{e)) = E e--^-r(^-) = E / e-^*P[r?. G dt, 7r(t) = i], \e\ < 1, 

is analytic inside the unit circle for all s, x > , then the function 
Q|(x)= i^-^>)Us-m) 

'^^ ' {i-\)f{s-k{e)) + \-B ' - ' M y> y) 
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is analytic on the open set \9\ < c(s). In this region it can be represented as a 
powers series 

Q^(x) = J] ^'^QKx), s,x>0, \e\<c{s). 

The coefficients of this expansion can be calculated by means of the inversion 
formula 

27rz 6"+^ (i - A)/(s - k{e)) + \-9 

We will call the sequence {Ql.{x)} ■> a; > 0, defined by the formula ([8]) the 
resolvent sequence of the process {Dx{t)}t>o- 

We now explain a probabilistic meaning of this sequence. Introduce a random 
sequence as follows: (see [26] ) 

n 

Xo(x) = 0, =7r(r/,)-(^, = + ^ X,', X„ = X„(0), 

where X = 7r(r/) — (5 G Z, {X, X^}, n G N is a sequence of i.i.d. random 
variables. We now define a right-continuous step process in the following way: 

The sample paths of the process are constant on the time intervals [Tin{x)-,rin+i{x)), 
n G Z"*" and there occur jumps at the instants r]n{x), n G N. These jumps have 
the same distribution as X = n^rj) — S, where n G {2,3,...}, and Xi{x) = 
'^(Vx) — S for n = 1. Here and in the sequel we will call the process {Sx{t)}t>o 
a semi-Markov random walk generated by the sequences {r]n{x)}, 
n G Let = sup„<j 5*0 (tt) be the supremum {S(j{t)}t>o- The generat- 
ing function of was found in [26]: 

B0St. = (l-/»)(^-c(.)) ^ ^ 

l-c{s) e-X-{l-X)f{s-k{9y 

where i^s is an exponential variable with parameter s > 0, independent from 
the process {Sx{t)}t>o ■ It follows from ([7]) and from the latter formula that for 
1^1 < c(s) 

Comparing the coefficients of 6^, /c G Z+ in both sides yields 

Qii^) = ^^jp: j;c(.)-'=-i5]E [e-^^7r(r?,.) = j] P[5+ = z - j]. 
J^^) i=0 j=0 
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Denote by n'^{r}x) € Z"*", s > a random variable given by its distribution: 

P [vr^fe) = A;] = J—B [e-^''^7^(7?,) = k] , k G Z+ . 
fxis) 

Then the previous equahty implies that 

which explains the probabilistic meaning of the resolvent sequence. Asymptoti- 
cally, one has that Qf,{x) ~ c(s)~'^ as A; — oo. 

Let Xq = {0,x}, X > 0, k £ Z+ and introduce upper one-boundary func- 
tionals of the process {Xt}t>o ■ 

t\x) = inf{t : D,{t) > k}, T\x) = D.,{t\x)) - k, r^^x) = r? + (r'=(x)) 

i.e. the instant of the first crossing of the level k by the process {-Dx(i)}t>0) the 
value of the overshoot across the upper level and the value of the linear component 
ry^(-) at the instant of the first crossing (the time since the last renewal). Denote 
*B^(x) = {r''(x) < oo}, 



f{x,dl,m,s) =E e-"^'(^);7?*^(x) £ dl,T''{x) = m,^''{x) 



m G N. 



We now determine the upper one-boundary functionals of the process {Dx{t)}t>o- 
Let /c G Z+ and = inf{t : 7r(t) > A;}, = ir^f^) - k he the first crossing 
time through the upper level k by the compound Poisson process {7r(t)}j>o and 
the value of the overshoot at this instant. Denote by 



Pk 



(t) = P[7r(t) = A;], ^^VfcW = Er(*) =e*'=(^\ \e\<l, 



k=0 



p'^idt) = P[f ^ G (it, T'^ = m] = ^ ^ Pi{t)'P[x = k-i + m]dt, m G N. 

i=0 

Lemma 3 ( [25j ). Let {Dx{t)}t>o be the difference of the compound Poisson 
process and the compound renewal process, 5 ~ 5'e(A), x > 0, A; G Z+ and 
{Q%{x)} k<=:'L+ 1 be the resolvent sequence of the process {Dx{t)}t>Q-, given by 

(i) the Laplace transforms of the joint distribution of {T^{x),r]''{x),T^{x)} sat- 
isfy the following equality 

f\x, dl, m, s) = e-^C--) > x}pT{d{l - x)) 

1 — r [x) 

k 

+ $1(0, dl, m) QUx) - e-^' [1 - F{1)] ^ Q|(x) p^_,(d/), (9) 

1=0 

where ^l{0,dl,m) = e-*'[l - F(/)] c(s)Vr(^^0; 
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(ii) the Laplace transform of the first crossing time through the upper level k 
by the process {Dx{t)}t>o are such that for all k G Z+, s,x >0 

A.,.)^Ee--^(^)^l--^fiM-A^(.), (10) 

where = E p^{s) [l - Q|_,(x)(l - X)~^] , pk{s) = s e-^' pk{t) dt; 

(iii) for E[x], E[r]] < oo and p < 1, t^{x) is a defective random variable and 

P[r'=(x) < oo] = 1 - (1 -/))(!- A)~^Qfc(x) < 1, A; G Z+ x > 0, 

where {Qk{x)}j^(zz+ , x > is the resolvent sequence of the process {Dx(t)}t>o, 
given by ^ for s = : 

^ . N 1 / d9 (I - X)fJ-k(9)) , , , ^ 

Qkix) = 6 -T-T ^ — , aG 0,c(0 ); 11 

' 27riJ\e\=aO''+^{l-X)f{-k{e)) + X-e' ^ ' V V ; 

if p > 1, then for all k G Z^, x > t^{x) is a proper random variable. 

Along with expression (jlip there exists another way to calculate Qk (x) , which 
is more applicable from practical point of view. We will now derive the recurrent 
formula for Qk{x). It follows from ([7]) for s,9 = that 

Qo{x) = (1 - A)(A + (1 - A)/o)^Vo(x), 

where for all A; G Z+ 

/•oo 

fk{x) = P [7T{r],) = k]= P[r?, G dt, n{t) = k], A = A(0). 

Jo 

Again, it follows from ([7|) for s = that 

(1 - x)U-k{e)) = (1 - X)f{-ki9))qe{x) + (A - e)qe{x). 

Comparing the coefficients of 6^, /c G N in both sides implies that 

k 

(1 - X)fk{x) = (1 - A) ^ Q^{x)fk-i + XQkix) - Qk-i{x). 

i=0 

Combining like terms yields 

k-l 

(A + (1 - A)/o) Qk{x) = (1 - X)fk{x) + Qk-i{x) - (1 - A) ^ Q^{x)fk-i. 

1=0 

The latter formula is a recurrent relation which allows to calculate successively 
the terms Qk{x) given the previous terms Qo{x), . . . ,Qk-i{x). For instance, 
given the expression for Qq{x) one finds that 

1 - A 
A + (1 - A)/o 

The knowledge of the one-boundary characteristics of the process allows us 
to solve the two-sided problems, which is the aim of the following section. 



Qi{x) 



^ , l-(l-A)/o 
+ A + (l-A)/o 
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2 Two-sided problems for the process {Dx{t)}t>o 

Let B G Z+ be fixed, k £ 0,B, r = B — k, Xq = {0, x}, x > 0, and introduce 
the random variable 



the first exit time from the interval [— r, A;] by the process {Dx{t)}t>o- This 
random variable takes values from a countable set {S,n, n G N} U { r]n{x) , n G N}, 
and it is a Markov time of the process {Xt}t>o {Ck are the instants of jumps of 
the process 7r(t). ) Note, that the exit from the interval can occur either through 
the upper boundary k, or through the lower boundary — r. In view of this remark 
introduce the events 

= {D^{x) > k}, i.e. the process {Dx{t)}t>o exits the interval [— r, A;] through 
the upper boundary k ; 

Sir = {Dx{x) < i-e. the process {Dx{t)}t>o exits the interval [— r, /c] 

through the lower boundary — r. Denote by 



T=(Z),.(x)-A:)% + (-D.(x)-r)l2,., L = ry + Wla^ + 0-%,, P[2t'=+2l,] = 1 



the value of the overshoot through the boundaries of the interval [— r, k] by the 
process {Dx{t)}t>o and the value of the linear component at the instant of the 
first exit (the time since the last renewal), where la = l2i(<^) is the indicator 
function of the event 21. Denote 



Theorem 1 ( |25] ). Let {Dx{t)}t>o be the difference of the compound Poisson 

process and the renewal process (2), 5 ~ 5'e(A), {<5|(a;)}fcgz+, he the resolvent 

sequence of the process given by Q^. =^ Qfc(O). Then 

(i) the Laplace transforms of the joint distribution of {x, L, T} satisfy the 
following equalities for all x, s > 0, m G N 




inf{t: aAt)^[-r,k]}''^'x 



V''{x,dl,m,s 



) = E e-'^;Ledl,T = m,^^ ,V;(x,m, 



s) = E[e-'^;T 



m, 2tr] . 





Bf^+^{0,dl,m,s 



(12) 



where the function f^{x,dl,m,s) is given by ^ 



Eg|+5 = j;(l-A)A'=-iQ^+5, 



Bf^+^{0, dl, m, s) = ^(1 - A)A'=- ^=+^(0, dl, m, s); 



fceN 
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(ii) for the Laplace transforms of the first exit time x from the interval by the 
process {Dx{t)}t>o the following formulae hold 

E[e--2l.]=|i^, (13) 



E 



where EAf+^{s) = E (1 " A)A^-iA^+^(s); 

feeN 

(iii) the probabilities of the exit from the interval through the upper and the lower 
boundary by the process {Dx{t)}t>o are given by 

P[2i,] = P[2t'] = 1 ^'^"'^ 



E Qs+B ' E Qs+B ' 

where the resolvent sequence of the process {Qk{x)}kez+ ^ x > is defined 
by m, Qk = Qk{0)- 

Denote by Ug ~ exp(s) an exponential random variable with parameter s > 
independent of the process Dx{t). For k G Z+, x > define D^{t) = supL'a;(-), 

D~[t) = inf Dx{-) the running maximum and minimum of the process. Our aim 

is to determine the joint distribution of {D^{t),Dx{t),D'^{t)}. In order to do 
this, we will require the joint distribution of {Dx{vs)-,D'^{i^s)}- 

Lemma 4. Let k e Z+ and E^{x,z,s) = E [z^^^"^'^; D+{us) < k] , \z\>l be 
the generating function of the joint distribution of {Dxii's), D^(i's)}- Then 

(i) the generating function E^[x,z,s) is such that 

k-l 

E^{x,z,s) = z'A'xis) + (1 - z)^z'Ai{s) + z'=Q|(x)E^/,(0, z), (14) 

i=0 

where 

^A/.(U, z) - l-c(s)/z' 

(ii) the joint distribution (B^[x,u, s) = P [Dx{i's) ^ u,D^{i's) < k] , n G — oo, A; 
satisfies the following equality 

^+{x,n,s)=A^x{s) + ^^j^^c{s)''-^QUx), Ali{s) = 0, u < 0; (15) 



(iii) under the condition (A) 



(A) p= (1 - A)^Er/Ex= 1, a'^ = fi Ex(x- 1) + E>rEr?- 



,2 



(l-A)(Er,)^ 



< OO, 
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the following limiting equality holds as B ^ oo, k > 



P[D,itB^) < [uB],Dt{tB^) < [kB]] 



1 



2k-u 



e '-'/^^''dv, 



(TV 2vrt J -u 

where [a] is the integer part of the number a, u < k. 

Proof. In view of the total probability law, homogeneity of the process Xt with 
respect to the first component, Markov property of r]i{x) we can write for the 
function (x, z, s), k G Z"^, x > the following equation 

E+{x,z,s) = s J e-"*P [t]^ >t]E ^z^W ; ^(t) < A:J dt+ 

+ / e-''^Y.^[7i^edu,7r{u)=v]z^Y.(^-X)X'-'z-^E+_^^^{0,z,s). (16) 

n 1 



v=0 



r=l 



Introduce the generating function Eg(x,z) = ^ 6''E'^{x,z,s), \6\ < 1. Multi- 

plying (jl6p by 9^ and summing over k G we derive the following equation 
for the function Eg(x,z) 



EUx,z) 



s 1 - Us - k{ze)) 



l-e s- k{z9) 

+ f^is - k{ze))^^ ~ ^'^ 



(17) 



E 



X-z 

Letting x = in the latter equation yields 
E^(0,.) '^-^ 



^/,(0,z)-E^(0,z) 



< 1, \z\ > 1. 



{l-X)f{s-k{ze)) + X-z6 



f{s - k{ze)) 



1 - A 



lE^/.(0,^) 



A 



S l-f{s-k{z0)) 

-9 s- k{z9) 



The function which enters the left-hand side of this equation is analytic in \6\ < 1. 
In view of Lemma [1] it has denominator of the right-hand side has a simple zero in 
9 = c{s)/z. Hence, the nominator of right-hand side should also have the simple 
zero. Letting 9 = c{s)/z in the nominator we find the function E^^^(0,z) 



s(l-A)-i l-c(s) 



s — k{c{s)) l — c{s)/z^ 



z > 1. 



Employing the definition of the resolvent ^ and substituting the expression for 
E^(0,z) into ([n]), we get 



E^(x, z) = Af{s) + 9 [-| Af (s) + Q^,(x)E^/,(0, z), 



(18) 



where 



Af(s) = £(z^)^^^(.) 



k=0 



s — k{z6 



l-z9 1 - A 



(x) 
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Using the definition of tlie resolvent ([8]) and comparing the coefficients of 6^, 
k € Z+ in both sides of (jlSp imphes that 



k-l 



E+{x, z, s) = z'^A'.is) + {l-z)Y^ z^Aiis) + z^Q|(x)E^/,(0, z), 



i=0 



i.e. the equahty (|14|) of the lemma. Comparing the coefficients of z*, i G —oo,k 
in both sides of the latter equality, we find 

P[D,ii^s) = i, D+{us) <k] = 

= Aiis) - Al-\s) + c{s)'-'QUx), i < k, 

where ^^(s) = 0, for i <0. The latter formula implies (jl5p . Denote e\{x,u^B) = 
P [D^itB"^) < [uB],D+{tB'^) < [kB]] . It is clear that 

lim / e~'^^e\{x,u,B)dt = - lim gf/™ (x, ftii?]). 

To proceed further, we need the following limiting equalities (see |26] ) 
cis/B"^) = 1 - B-^y/%s/(j + o(fi-^), 



lim B-^Q'l^. (x) = ^ J- ^ ^ = lim S 



B-X^^ ^my-^^- ^^Er? "b^^^" 

lim A]^^\s/B'^) = 1 - ch f kV2^/a) = lim Ai^^''^\s/B^). (19) 



In view of these equalities and of the formula (jlSp we derive 

/ e-^*4(x, n, dt = s-^I{„<o} ( e"^/'^/2 - e-(2fc-«)v^/'^/2 
s"'l{«6[o,A.]} (l - e-^^/72 - e-(2'=-«)^/72) , u<k. 



lim 

B^oo 



Denote by w^t>o} the symmetric Wiener process with the dispersion a and by 
= inf{t : wt > a} the first passage time of the level a E M+. The Levy formula 
F [t < t] = 2P [wt > a] implies for the Laplace transforms that 



/•oo 

2 / e~'*P [wt > a] dt. 
Jo 



Employing the latter formula to invert the Laplace transforms in the previous 
equality, we derive the second limiting formula of the theorem. □ 



Let k,r G Z+, u € —r,k and denote by 

e*,fc(x,u) =P[-r< D-{t), {t) < n, Z)+(t) < A;] = P [Z),(t) < u, (r) > t] 

u) = e^/j^ix, u) = s e-'^el^^ix, u) dt 
Jo 

the joint distribution of {D^{t),Dx{t),D^{t)} and its Laplace transform. 
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Theorem 2. Let ~ exp(s) be an exponential random variable independent 
of the process Dx{t), B = r + k. Then 

(i) the joint distribution of {D~{us),Dx{i^s),D^{i/s)} is such that 



ek{x,u)=A:{s) 



EQ 



S+r+u I 



u G — r, k, 



(20) 



S+B 



where EA^+^+"(s) = (1 - A) £ A^-M^+^+"(s); 

i=l 



(ii) under the condition [A) and r € (0, 1), k = 1 — r the joint distribution 



■[vB] 



[fc_B](^' t^-^]) weakly converges as B ^ oo to the joint distribution 



-r < inf w^, Wt < u, supWy < k 



, u G [-r, k] 



of the infimum, the supremum and the value of the symmetric Wiener pro- 
cess with the dispersion a. In addition, the following limiting equality holds 

-4(7rn<T)2 



e 2\ 



n 



sin(r7rn) sin 



2 r + u 



■ nTT 



(21) 



Proof. The total probability law, homogeneity of the process Xt with respect to 
the first component, Markov property of Xri^) k,r G Z"*", x > imply 

the following equation for l^j > 1 



E £»+(/.,) < A:] = E (x) > 

oo 

+ J]]E [e-^>^^(^);r = i,2l,] z-('^+^)e [z^^^^^^; D+(i/,) < i + B 



(22) 



where B = k-\-r. This equation for the case of a spectrally one-sided Levy process 
was derived in [22] , and for the general Levy process in [29] . Let us briefly explain 
the equation (|22p . The increments of the process Dx{t) on the interval [0,1/^] 
without the intersection of the level k (the left-hand side) can be realized either 
on the sample paths of the process which do not cross the negative level — r (the 
first term of the right-hand side) or on the sample paths which do cross the level 
—r and then the further evolution of the process is nothing but its probabilistic 
copy on [0, Vs] (the second term). In view of ([22]) and (fT2]) . (fni) we find for the 
function E^f^{x,z) = E [z^^^*"'); (a:) > i^s] that 

E:.^,{x,z) = E^{x,z,s) - :^|^(1 - A)^A-iz-(^+^)i?+^(0,.,.) = 



k-1 



,''A'x{s) + {l-z)^z^Al{s)+ 



i=0 



+ z 



, Qlix) (1 - X)A^{s, A) - (1 - z){X/z)^+^A^{s, z) 



Q%{\) 



l-X/z 



(23) 
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B 



where = ^ zM^(s), Q%{X) = ^ A^Qf(O). The formula ^ yields 

j=0 i=B+l 

P [xf (x) > = A',{s) + ^M^ A^{s,X). 

It is not difficult to derive the following equality 
k 

J2 z''KA^,u) = ^ {E^rA^,z) - z'^+ip [xf (x) > i 



The right-hand side of (j23p implies that 



fc fc / N B B—i 

j=0 



M=-r M=0 QbW i=0 



Comparing the coefficients of z", u G —r,k, we find 
Since 



Qb\^) 



^QUb = (1 - A)A-«-ig^(A), AMj,(s) = 0, 



i=0 



one can see that the previous equality is the formula (I20p . Let us verify (I2ip . It 
is clear that 



rB],[fcB](2^' = ^[/^,[fcij](2;J^^^]), /i:G(0, 1) r=l-k, 
where the function is determined by (|20l) . Thus, 





1 -ch ( —V2^ 



1 sh/cy^/c 



s shV2s/o- 



ch( !:±^^/2^) -1 



, (24) 



where u'^ = max(0, u). In order to compute this limit we used the formulae (119 
Note, that the inversion of the Laplace transform in the right-hand side of ([2 
this equality was found in [22j and resulted into the following formula (a > 0) 



-r < inf Wy, Wt < u, supif^ < k 



1 



r>Q+ioo 

-/ e-e*{s)ds 

'' J a— too 



E 



g-|{7rna)2 



sm rvrn sm 



2 f r + u 



IT ^ — ' n 

neN 



■ nvr 



n G [— r, /c]. 



Therefore, we established the weak convergence of the joint distribution e* (n, B) 
as -B — ?■ CO to the corresponding distribution of the Wiener process and also 
verified the formula (1211). □ 
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3 Reflections from the boundary 



Denote by D'^[t) = r + Dx{t), t > the process starting from r E Z when 
Vxi^) = X > 0. Let B G Z"*" and for all i > we define right-continuous 
processes reflected at the boundary B as follows 



D^^ {x,t)=Dl{t)- max <0, sup Dl{-) - B} e-oo,B, r€-oo,B. (25) 

i [o,t] J 

The first reflection from the upper boundary B of the process D^. (x, t) takes 
place at r^~^(x). Then the process stays at the boundary for some random time 
rji, where / = r]^{T^~^{x)). At the instant t = t^~^{x) + rji the process is 
reflected to a random state B — 6. In the sequel the evolution of the process 
) is a probabilistic copy of its evolution on [0, '''{x) + rji). It is worth 
noticing that reflections from the boundaries reflected by infimum (supremum) 
were introduced by Levy for a standard Wiener process. Applying the symmetry 
principle and the mirror reflection principle Levy determined the distributions of 
the boundary functionals of the reflected standard Wiener process. We will show 
that these distributions are the weak limit distributions for the reflected process 
after an appropriate scaling of time and space. 

3.1 Passage of the lower boundary 



We now define the boundary functionals for process (|25]) . For r G 0, S denote 

rf (x) = inf{t : Df{x, t) < 0} = r, rf (x) = -Dfir) = T, r G [0, B] 

the first crossing time of the lower level by the process D^. (x, t) and the value 
of the overshoot at this instant. Note, that these boundary functionals were stud- 
ied in [3n| for the reflected Levy processes generated by infimum (supremum). It 
is worth noticing that in this article the asymptotic expansions for the distribu- 
tions of the characteristics of the process were determined for the reflected Levy 
processes obeying the two-boundary Cramer's conditions. 

The reflected spectrally one-sided Levy processes generated by the infimum 
(supremum) of the process were considered in [3], [41]. An interesting applica- 
tion in queueing theory for the spectrally one-sided Levy process reflected by its 
infimum was given in [Bj. 

Theorem 3. Let {D^. {x,t)}t>o be the reflected processes defined by [25\) . B G 
Z+, rG(VB, 

y^(x, dl, m,s)=B \e-'^; L€dl,T = m, 2t^'j , K(x, m, s) = E [e""^; T = m, 21^] 

the Laplace transforms of the joint distribution of {Xr{x), L, T} of the process 
{D^{t)}t>o Then 
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(i) i/ (5 G N (an arbitrarily distributed non-negative variable), then the Laplace 
transform of the joint distribution of {r, T} is such that for (m € N) 



vi(x 



,m) = E e-'^-^''^;T = m^ =Vr{x,m,s) 



1 - ^(0) 



P[5 = m + B] + Y, = i]VB-i{0, m, s) 

i=l 

where V^{x, dl, s) = Ylm=i ^'^{x, dl, m, s), 

a\x)= V\x,dl,s)fi{s), A{x) = y^V[5 = k]a\x)- 
■^0 k=i 



(26) 



(ii) if 5 ^ 5e(A), then the following equalities hold 
fl{s) + {l-f{s))SI,_,^_,{x) 



vf.{x, m) 



/» + (l-/(s))E5|^^_i 



(l-A)A™-\ re 0,5, (27) 



where Sl{x) = Y.toQt{^)^ E 5|^^_i_= (1 - A) E.=i A'-^5f_i+5(0); the 
random variable t^{x) is proper (P [r^(x) < oo] = 1) and 

Erf (x) = Er/^ -Brj + Br] [BSs+b-i - SB-r-i{x)] < oo, (28) 
where Sk{x) = S'^i^x), ES'^+b = E5^_^^; 
(iii) under the conditions (A) the following equality is valid 

lim Ee [■'^1^ " = — r G (0, 1), k = 1 — r. 



ch (\/2i/a) 



Proof. Let us verify the formula (I26p . It follows from the definition of the process 
Dj. {x,t) (p5]) . the total probability law and the Markov property of Vn{x) 
that the following system of the linear integral equations holds 

POO 

v^{x,m) = Vr{x,m, s) + / V^{x,dl,s)v%{l,m), k = B — r, 



m] 



v%{x,m) = fl{s)P[5 = m + B] + /,(s) J^P[5 = r]v%_,{0, 



r=l 



This system is similar to a system of linear equations with two unknowns and 
can be solved analogously. Substituting the expression for v'^^{x,m) from the 
second equation into the first one, we find that 



l{x, m) = Vr{x, m, s) + a''{x)P[6 = m -\- B] + a^{x) ^ P[(5 = r]w^_^(0, 



m). 



r=l 
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Letting a; = in the latter equation after calculations yields 



J2 ^[^ = r]v%_riO, m) = -P[5 = m + B] 
B 

P[6 = m + B] + Y^ P[6 = k]VB-k{0, du, s) 



+ 



k=l 



il-AiO))-\ 



Inserting the right-hand side of this quality in the previous one, we get (j26p . 
In case when 6 ~ ge{X) the formula (I26p takes a more simple form. The first 



formula of (fT2]) and (f26]l imply that (x) ~ ge{X) for any r G 0,B. Summing 
over m G N both sides of (j26p . we find for the function v^{x) = Ee"'*'^'' that 



v^r{x) = Vr{x,s) + 



a^-"(a:) 



B 



(29) 



Now we calculate a*^(x), ^(0) in case when 5 ~ ge(A). Employing the formulae 
(|12p and performing the necessary calculations we find that 



aHx) = Us) + {l-f{s))Si_,{x) 
1 - A 



5+B 



1-^(0) 



B 

A^ + (l-A)^ Wb_,(0,s) 



l-Qo(O)] [/(s) + (1 - /(s))E5|+5_i 
1 - A 



i=l 



XEQ 



[i-Qo(o)] 



5+B 



Substituting the right-hand sides of these equalities into (j29|) and taking into 



7^B , 



account that (x) ~ 5'e(A), we derive (f27|l . 

The limiting formulae (|19p were obtained in [26j. Similarly (see [26]) we derive 
for all A; > that 

lim S'^l^' (x) = (ch (kV2^/a) - l) = lim B'^ ^S'IbLs' (30) 



sE?7 



'[fc_B]+<5- 



Letting B — t- oo, we have fx{s/B'^) = 1 — Eijxs/B'^ + o{s/B'^), which implies 
for r G (0, 1) that 

lim Ee-^^?s](")/^' = 1 + - 1)) ^ ch {k^/a) ^ ^ ^ _ ^_ 



1 (ch (^/2i/cj) - 1)) ch(^/2^/o-) 



The equality (I28p follows from the following ET^(a 



□ 



3.2 Increments of the process reflected in its supremum 



Define DQ{x,t) = ^^^(t)— max< 0, supL'a;(-) — /c > G —oo,k, the process refiected 

i [o,i] J 
from the upper boundary k G Z"*" generated by its supremum. 
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Theorem 4. Let {DQ{x,t)^t>o be the process reflected from the upper bound- 

, M S —oo,k be the distribution of its 



ary and p\{x,u) 



Dq{x,Vs) < u 



increments on the exponential interval [0,z^s]. Then 



(i) for all k G Z+ k) = 1; and for u G -co, /c — 1, x > the following 

equality holds 



pI{x,u)= Al{s) + c{s) 



k—u—l 



fAs) 



(31) 



l-f{s) 

where A^{s) = 0, for u<0, F{s) = s{l - c(s))/(l - A)(s - k{c{s))); 
(ii) under the conditions (A) for k > 0, u < k the following relation is valid 



lim P 



D^o^\x,tB^) < [uB] 



1 



1 



(T\/2vrt Ju 



2k-u 



e "'^'^"^'dv; (32) 



(iii) if p > 1, then the ergodic distribution Pk{u) 
ists, 

Ex 1 - c 



lim P 



D^{x, t) <u 



ex- 



Pk{u) 



Jt — U—l 



p 1-Ec* 



, u E — cx), k — 1, c = lim c{s) G (A, 1) 



Proof. Define the generating function distribution of the process 

k 



pI{x,z) = Ez^o(^.^.) = ^z'P Dq{x,u, 



It is obvious that P^(x,l) 



\^Dq{x,1^s 



< k 



\z\ >l,k£ Z+. 



1. In accordance with the 



total probability law and the definition of the process DQ{x,t) we can write 



pUx,z) = E+{x,z,s) + z^ 



/|(d/)(l-/z(s)) 



+ z' 



poo 

/ f^{dl)Ms){l - A)^ A-1z-P:(0,z), 

■^0 ^=l 



(33) 



where the function E'^{x,z,s) = E [e t'^'(x) > i^s] is given by and 

the function f^{dl) = E e~^'^''^^^;r]'^{x) G dl is determined by ([9]). This equa- 
tion means the following. The sample paths on which the increments of the 
process Dx{t) occur can be decomposed into three types: 1) the sample paths 
which do not intersect the upper boundary k (the first term of the right-hand 
side); 2) the sample paths which do intersect the upper boundary and stay there 
(the second term); 3) the sample paths cross the upper boundary and then they 
are refiected (the third term). After some calculation which we skip, one can see 
that the formula ([9]) implies that 

f^idUis) = Us) + (1 - f{s))Si{x) - -— ^ Ql{x). 
J- - c{s} 
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Letting x = in ()33p and taking into account the latter equality, we derive 

oo 

■ , ■ o ^ ( .s ) \ / 7. — \ 

z > 1. 

i=l 



(1 - A) V X'-'z-^PtiO, z) = , 



Substituting the right-hand of this equality and the expression (I14p for E^{x, z, s) 
into das]), we find that (|z| > 1) 



Pl{x, z) = z'' + {I - z)y z'Aiis) + 2^(3) ^% ^] ( -MfL 

t'o ^l-cis)/z \l-f{s) 

Comparing the coefficients of z*, i £ {k,k — 1, . . .}, we get 



P 
P 



DQ{x,iys) = k 



l-At\s)-F{s) 

D',{x,u,) = i]=Ai{s)-Ai7\s)+ 

Us) 



+ F{s)cisy-'-'il-cis)) 



1 - fis) 



i < k. 



One can see that the second formula implies the equality (|3ip of the theorem. 
Let us verify (|32p . It follows from the first formula of ()19p that 



F(s/B^) = -^Brj + oiB~^), lim c(s/52)[B(fc-«)]-i ^ g_(fc-«)v^/a^ 



Denote p^f.{x,u,B) = P D^Q^\x,tB'^) < [uB] , k > 0. Then ([E]), 
lim / e ^ pdx,u,B) dt = - lim PrtRi (a:, [ni?]) = 

B-s>oo Jq S _B-!>oo Lf^-DJ 

^ (l - ch [u+V2's/a) + e-('=-«)v^/'^ch (fcV2i/cJ^ 



s 

_ 1 
s 



^{'u<0} 



^uV2s/a j2 _|_ g-(2fc-M)V2s/o- _)_ 



{«e[o,fc]} 

where = max{0, u}. Employing the formula ^ g-av^/o- _ 2 J^°^ e~'^*P \wt > a] dt, 
to invert the Laplace transform, we derive the limiting equality of the theorem. 

For /O > 1 the mathematical expectation of t^{x) is finite. It follows from 
(fTOl) that 



Et*^(x) 



(1 - \)k{c) 



Pi 



i=0 



1 



Qk-i{x) 

1-A 



< 00, 



where pi = lim s ^Pi{s) = P [7:{t) = i]dt < 00. Moreover, the process 



s-s>0 

DQ{x,t) is of regenerative type [33j. The instants if the passages of the upper 
boundary are the regeneration times. Hence, [33] there exists the ergodic distri- 
bution of the process Pk{u) = hm P D^{x,t) < n . To determine this distribu- 

f— ^00 L 

tion, it suffices to apply to (j3T]) the Tauberian theorem: pk{u) = \mipi{x, u). □ 

s— >0 
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Let {DQ{x,t)}t>o be the process reflected from the upper boundary. Define 

for r,k & Z+ 



Tr,k{x) = inf{t : DQ{x,t) < -r} =^ r, 



Tr,k{x) 



the first exit time from the interval [— r, /c] by the process DQ{x,t) and the value 
of the overshoot through the lower boundary — r. Since Xt is homogeneous with 
respect to the first component, then {Tr,kix),Tr^kix)} are identically distributed 
as {t^{x),T^ (x)}, B = k + r and their joint distribution is determined by ([27j) . 
^ 

Theorem 5. Let {Dq{x, t)}t>o be the process reflected from the upper boundary, 



Pr.k{x,u) 



u ^ —r,k be the distribution of the 



Dq{x,Us) < u; Tr^k{x) > V. 
increments of the process on the interval [0, z^^] on the event {T^^k{x) > Vs}. 

(i) the distribution of the increments is such that for all r,k £ Z^, 
fl{s) + {l-f{s))Sl_,{x) 



Pr,k{^^ /c) 
Pr,kix,u) 



f{s) + {l-f{s))BSI^^_, 
. _ fl{s) + {l-f{s))Sl_,{x) 
"^'^ /» + (l_/(s))E5|_,^_, 



B = k + r, 



(34) 



e4+"+"(s), ue-r,k-l; 



(ii) under the conditions (A) the following limiting equality holds 



lim P 

S— >-oo 



[do''\x,tB^) < [UB]- T^rBUkB]{x) > tB' 



dcf 



- E 



e 2^ 



vr 



n + 



sm 



pit) -. 
r + u 



(35) 



n H — vr 
2 ' 



where r G (0, 1), k = 1 — r, u £ [—r, k]. 

Proof. In accordance with the total probability law, homogeneity of the process 
Xt with respect to the first component, Markov property of Tr^ii.{x) and the 
properties of the exponential variable Vg we can write 



Pk{x,u) = p^^k{x,u) +vf.{x){l - X)'^X' V?+B(0,M + r + i), ue-oo,k, 



i=l 



where the function pl.{x,u) 



Dq{x,Us) < u 



is determined in Theorem HI 



-r^k, 



This equation means that increments of the process Dq{x, Ug) are realized either 
on the sample paths which do not exit the interval [— r, fc], or on the sample 
paths which do exit the interval and the further evaluation of the process is 
its probabilistic replica on [0, i^^]. Substituting the expression for the function 
pI{x,u) into (f3T]) after necessary calculations we derive 
For r G (0, 1), k = 1 — r, u G [— r, k] denote 

pi^kix,u,B) = P M^''\x,tB^) < [uB]; T[,BUkB]{x) > tB' 
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Employing the third formula of ()19p , the limiting equality of Theorem [3l we find 
1 /r2 f°° 

- hm pj^/l, {x, [uB]) = hm / e-'%^k{x, u, B) dt = (36) 

l-di{u+V2^/a) 1 ch{kV2^/a) / r— \ dcf . 

= + - , , , (ch (n + r V2s/a) - 1] = p s), 

s s ch(v2s/cr) V / 

where = max{0,u}. When u G [— 0] we derive from this formula that 



- v.; , sh^ V2^/a , ue[-r,Q]. 

s ch(V2s/cr) V 2 y 



It is clear that s = is not a singular point (pole or point of branching) of the 
function p*{s). In the semi-plane 5f?(s) < this function has simple poles in 

1 2 2 / ^+ 
Sr, = crvr nH — , nG Z , 

2 V 2y ' 

and it is analytic in the whole plane apart from these points. Hence, for a > 
p{t) = — / eV(s) = Res,=,„/(5). 

J a— too 

Calculating the residues of the function p*{s) in s^, we obtain the right-hand 
side of the formula (j35p for G [— r, 0]. On can see that the first term in the 
right-hand side of (I36p is analytic in the whole plane for u G (0, A;]. Applying the 
inversion formula, we find that the contour integral of this term is equal to zero. 
The second term of ()36p is the same also for u G [— r, 0]. Thus, the formula (|35p 
holds for u G [— r, k]. □ 

4 Applications forM^|G^|l|B system 



Let B G 1r , r G 0, + 1, x > 0. Introduce the two-component Markov process 



Yr,x{t) = {4,x(t),??r,x(i)} G 0,5 + 1 X M+, Yr,M = {r,x) 
by means of the following recurrent equations 

^"^^"Uo,o(t-f/^-^Hx))/ t>f--Hx), ' + 



(0,0), < t < /I ~ exp(/i). 



lo,o(t) = <( Yr,o{t - /i) : P[x = r], r = 1, 5, t > /i, (38) 
YB+i,o(t-/i) :P[x>S + l], t>/i. 



where f/^+^(x) = inf{t : Df~^^{x,t) < 1}, r = 1,5 + 1. 
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Remark 1. Since the process Xt is homogeneous with respect to the first com- 
ponent, then the random variable f^~^^{x) is identically distributed as 



('^7|j and, hence. 



,-sf^+''{x)~B+l 



e 1,5 + 1. 



The process ^,a;(i){t>o} serves as a mathematical model of the functioning 
of M^|G^|1|B, (6 ~ 5'e(A)) system, which has the following properties 

(i) The customers arrive in groups (batch arrivals) according to the Poisson 
process with intensity > 0. The number of the customers in every group 
is represented by the random variable x G N. 

(ii) The system has a finite waiting room (buffer) whose size is equal to B + 
1 < OO. Suppose that upon the arrival of a new claim of size x it finds 
r € 0,B + 1 occupied space in the waiting room. Then min{A;,x} joins 
the queue, and loss of size max{0, >c — k} occurs, where k = B + 1 — r is 
the size of empty space in the waiting room; 



(iii) 



The duration of service completion is arbitrary distributed as 77 > 0. Sup- 
pose, that at time t the service cycle is accomplished. Then the occupied 
space in the buffer is reduced by min{r, 5}, where r € 1, .B + 1 is the value 
of occupied space in the waiting room at time t — 0. If at the instant of 
the service completion r — min{r, 6} > 0, then a new service cycle starts. 
If at the instant of the service completion r — min{r, 5} = 0, then the new 
service cycle starts upon arrival of a new claim (after exponential time with 
parameter ^ > ) . 



For all t > the event {Y^^xit) = {i, y)}, i £ 1, B + 1, y > means that at 
time t there are i occupied places in the waiting room, and y stands for time 
elapsed since the beginning of the service cycle. Here (r, x) is an initial state of 
the system. 

The event {Yr^x{t) = (0,0)} means that at time t the waiting room is empty 
and the server is idle. The system stays in the (0, 0) state for an exponential 
period of time (with parameter /i.) 

Therefore, dr^x{t) is the number of the customers in the waiting room at time 
t. If dr^x{t) > 0, then rjr^x{t) is the time elapsed since the last start of the service 
cycle up to time t. If d^^xit) = 0, then P[r7,.^i.(t) = 0] = 1. 



4.1 Busy period of the system 

Assume that at time to = system is in the state (r, x), where r E l,i? + 1 is 
the number of the customers in the waiting room, and x > is the duration of 
the current service cycle. Introduce the random variable 

hr{x) = inf{t : dr^xif) = 0} 

i.e. the instant at which the system for the first time becomes empty. Thus, the 
interval [0,6r.(x)] is a busy period of (r, x) type. 
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Corollary 1. Let hl{x) = E [e ^'^^^^'> ]br{x) < oo] he the Laplace transform of 
the busy period of type {r,x). Then the following relation holds 

fl{s) + {l-f{s))SI,_,{x) 



b'ix) = ^^ ^ , ' " ^-''^ ' , x>0, rG 1,5 + 1, (39) 

f{s) + {l-f{s))BSI^^_,' - ' 

where 

k oo 

5|(x) = J2Qt{x), E5|+5_i = (1 - A) 5] V-i5ti+5(0). 

i=0 i=l 

Observe, that the random variable br{x) is proper (P [br{x) < oo] = 1) and 

E6,(x) = Er/^ -Bt] + Bt] [BSs+b-i - Sb-Ax)] < oo, (40) 
where Sk{x) = Slix), BSs+b = 'ES^^^. 

These formulae follow straightforwardly from Theorem [3] and Remark [TJ 

4.2 Time of the first loss of a customer 

Suppose that the system starts functioning from the state (r, x) and denote by 
lr{x) the time of the first loss of the customer (a group of customers). 

Corollary 2. Let lf.{x) = E [e"*''"^^'); /^(x) < oo] be the Laplace transform of 
lr{x). Then the following relation is valid 



m+B - jfji Qi^) 

m+B - 7^, Qi^) 



BAf+^(s) - -^A(s) - 



where rSl,i? + l, k = B + 1 — r, 

B+l B+l 



i=l i=l 



Note, that the random variables /o(0), lr{x) are proper, and they have finite 
mathematical expectations. 

Proof. Let r G 1,S + 1, x > 0. Denote by x?^^{x) = inf{t : r + D^t) ^ [1,5 + 
1]} the first exit time from the interval [1,5 + 1] by the process r^D^it). Since 
the process Xt is homogeneous with respect to the first component, then the 
random variable y^^^{x) is identically distributed as x^_i(a;) and its Laplace 
transforms are determined by the formulae of Theorem [TJ In accordance with 



23 



the definition of the process we can write the following system of the 

equations for the functions lf.{x), Iq{0) 

l'^{x) = s) + Vr-i{x, s) /g(0), r e 1,5 + 1, X > 0, 

B+l 

m = -^as^i + E (42) 

S ~\~ Ll S ~r M . 

1=1 

where Oj = P [x = i] , Oj = P [x > , and 

y^_i(x,s) = — — , (43) 

V^^'-^ix,s) = 1 - A^-^'-^is) - (l - E , 

Substituting the right-hand side of the first equation of ()42p for x = into the 
second one, we get 

B+l B+l 
i=l 1=1 

The latter equation yields 

/ B+l \ / B+l \ 

m = J^UB+i + Y^^a,V''+'-'{0,s)\ h --^^a.F,_i(0,s)j . 



Taking into account the first formula of (I42p and (I43p . we derive the equalities 
[TJ of the corollary. □ 



4.3 Number of the customers in the system 

Let ~ exp(s) be the exponential random variable with parameter s > 0. 
Introduce the transient probabilities of the process dr^xit){t>o} ■ 



q^^^{u) = P[dr^xii^s) < u], Qo^Qiu) = P[dofi{us) < u], r,u£l,B + l, 
<,(0) = P[dr,M = 0], g^,o(0) = F[doA^s) = 0]. 
Denote b{s) = as &|j+i(0) + Ef=i 6|(0). 

Theorem 6. The distribution of the number of the customers at time Vg is 



such that 



lU^) = E<"-i(s) + ^^"^''^j^ . ql,{B + 1) = 1 - " 

s + /i — /io(s) s + /X — ^o(sj 

qUu) = Al-^{s) + bl{x) ^"^"'^L . + 1) = 1 '^^^''^ 



s + fj. — fib{s) ' s + n — iJ,b{s) 

s bf.{x) 



?o,o(0) = — 7-, <.(0) 

s + /i — /io(sj s + fj, — n b[s) 

where 



Cu{s, X) = ^-s + X{s + fi) [a-{s) - BAf+^{s)) . 
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Corollary 3. Let vrj = lim P[(i(.)(t) = i], i G 0,5 + 1 he the stationary 
distribution of the number of customers in the system M^IC^jllB. Then 



-, -1 



ttb+i = 1-M^ + Cb{X)), = 7ro(C,(A)-Q_i(A)), i£l,B, 

where pi = lim s~-'^/5j(s) = P [7r(t) = i] < oo, Co(A) = 0, 



C.(A) = t%-1+Am (a;^ - E^^+") , yl^ = lim -^I^(s) = f^/^, 
i — A \ / s— >0 s ^ — ' 

i=0 



Qu—i 

1 - A 



Proof. By pf.^{u) = P [(ir,x(i^s) < ftr(a;) > J^s] , r, n G 1, S + 1 denote the tran- 
sient probability of the process dr^xi^^s) on the event {^^(x) > I's}- Taking into 
account the homogeneity of the process Xt with respect to the first component, 
the definition of the process dr^x{t) and the formulae p4p of Theorem [5l we 
derive 



p'^^^iB + 1) = 1 - Pl^^iu) = Al-''{s) - hl{x)'EA, 



> 



n G 1,5, 



where the function is determined by (I39p . 

In accordance with the definition of the process l^r.^^lO we can write the 
following equations for the functions qf.^{u), Qooiu) for u G 1,5 



s + fl 



B 



aB QB+ifiiu) + ^a'iqlo{u) 



1=1 



(44) 



Substituting the right-hand side of the second equation into the first one, we get 



(u) =plJu) + hl{x] 



s + /X 



]{s,u), 



where q{s,u) = ob q^+i oi'^^-) + J2f=i o-i it oi'^) ■ Letting x = in the latter 
equality implies that 

aBqh+ifiiu) = dBP%+i,o{u) + aB 6|j+i(0)^— g(s, u), 

S ] fJ/ 



B 



B 



B 



^aiqloiu) = "^a^plQ{u) + bt{0)——q{s, 

i=l 1=1 i=l ^ 

Adding these equalities, we obtain the function q{s, u) 



u . 



j{s,u) =p{s,u) 1 



S + /U 



-1 



h{s)\ , nGl,5, 
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where b{s) = 0.3 6|+i(0) + E^ti a* (0) 



B 



p{s,u) = aBfB+^,,{u) + ^aipl,{u) = ^aaV(.^) " Hs)^4+--\s) 



4=1 



1 = 1 



Substituting the expression for the function q{s,u) into ()44p . we find that 



' (n) = Ar'-(.) + 6^?(x) 



Cu(g, A) 
s + ^ — fJ-b{s) ' 



s + fj, — fi b{s) 



(45) 



where 



1-A 



To derive the latter equahties, we used the fohowing relation 

u 

^i^ai A^-\s) -{s + fi)A^{s) = Qi-s, ue T;^, 

i=l 

which follows from the definition of the function A^{s). If u = B + 1, then 
Prxi^ + 1) = 1 — bf.{x), p{s, i? + 1) = 1 — b{s) and the following formulae are 
valid 

qloiB + l) = l 



■JB + 1) = 1 



'fix) 



(46) 



s + fi — fib{s) ' s + 11 — /u6(s) 

Taking into account the definition of the process Yr^xit), we can write the fol 
lowing equations for the functions q^^^{0), 9o,o(0) 

g,%(0) = bf{x)q'o,o{0), 



s u 

+ 



S + fl S + IJ. 

Solving this system yields 
9o,o(0) 



«i?9B+i,o(0) + X]"^<o(0) 



i=l 



Qr,M 



s bf.{x) 



(47) 



s + — fj,b{s) ' s + iJ, — fib{s) 

Observe that lim74^(s) = limEAo+"(s) = 0, lim6^(x) = 1. It follows from 
^- (|47p and properties of the Laplace transforms that 



limg; (u) = linigo,o(^i) = q{u) = lim P[(i(.)(t) < u], uel,B + l, 
lim (7,%(0) = lim g^^o(O) = ^(0) = lim P[d(.)(t) = 0]. 



The formulae (147^ imply that 



g(0) = lim 



s + ^ — [ibis) 
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In view of (gS]), ([M]) we find 



q{B + l) = l-q{0), qiu) = qiO)CuiX), uGl,B, 
where pi = lim s^^pi{s) = P [7r(t) = i]dt < oo, 



i=0 



^ Qu- 



1 - A 
□ 



4.4 M^|G|1|B system 

Let us stress the fohowing fact. If we set the parameter A = in the geometrical 
distribution P[(5 = n] = (1 — A)A"~^, n S N, A G [0, 1) of the random variable 
(5, then Y'[5 = !] = !. In other words it means that the process {Dx{t)}t>Q has 
unit negative jumps at the times instants {rin{x)}n&i and 5Ar^(() = Nx{t). Then, 
it follows from ([2]) that 

D:,{t) = Tr{t) - N^{t) G Z, t > 0. (48) 

We will call this process a difference of the compound Poisson process and a simple 
renewal process. Setting the parameter A = in the statements of Lemma [1] 
leads to the following result. 

Lemma 5. For s > the equation 6 = f{s — k{6)) 

has a unique solution c{s) inside the circle \6\ < 1. This solution is positive, 
c(s) e (0,1). // E[xl,Ey < oo, p = ;uE[xlE[?7l, then for p > I, lim c(s) = 

c G (0, 1); and for p < 1, lim c(s) = 1. 

The statements of Lemma's [2H1] and Theorem's [iHS] can be reformulated in a 
similar way. Letting A = in the defining formula ([8]) for all s,x > we get 

Q|(x) = -L/ 1 Us-m) aG(0,c(s)) (49) 

^'^^ 2m e^+^ f(s - k{9)) -9 

the resolvent sequence of the process {Dx{t)}t>o-, which is given by (PS]) . This 
resolvent sequence has been introduced in [23]. Setting A = in ([9|), ([TO]) , we 
obtain 

f^Adhm,s) = e-^('-^) > ^}pi:{d{l - x)) 

k 

+ %{dl, m) Qlix) - e-«' [1 - F{1)] ^ Q|(x) p^_,{dl), 

1=0 

Ee— '^W = 1 Qlix) - A^is) (50) 

S — K[C[S)) 
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the Laplace transforms of the upper one-boundary functionals of the process 
{D^{t)}t>o dH, where p,{s) = s e-«*P[^(i) = i] dt, 

k 

Al{s) = Y.h{s) [1 - QUix)] , nidl^m) = e-^\l - F{1)] c{sfPk{dl). 
i=o fcez+ 

We have introduced the auxihary functions and the resolvent sequence of the 
process (|48p . therefore we can state the following result. 



Corollary 4. Let {Dx{t)}t>o be the difference of the compound Poisson process 
and the renewal process (4^ , {Qk(.x)}k&z+^ x > be the resolvent sequence of 
the process given by ( f^PP , Qf, =^ <3|(0). Then 

(i) the Laplace transforms V^{m,s), Vj!^{dl,m, s) of the joint distribution of 
{x, L, T} satisfy the following equalities for all x, s > 0, m G N 

Vr{x, i, s) = Sa, VHx, dl, i, s) = f\x, dl, i, s) - ^|M/B+1(o, dl, i, s), 

where 6ij is the Kronecker symbol and f^{x,dl,m,s) is given by / fJQp ; 

(ii) for the Laplace transforms of the first exit time x from the interval by the 
process {Dx{t)}t>o the formulae hold 

E [e-x.2i,] = E [e-^>^; 21^=1 = 1 - A^s) + ^(1 - A^+^s)); 

(iii) the exit probabilities from the interval by the process {Dx{t)}t>o satisfy the 
equalities 

^B+1 UiB+1 

where the resolvent sequence of the process {Qk{x)}k&+ ■• x > Qk *== 
(5fc(0) is given by for s = 0. 

In order to prove the corollary, one has to put A = in the statements of 
Theorem [TJ The results obtained in Corollary U] can be applied for studying 
the queueing systems M^|G|1|B (P[5 = !] = !) with finite waiting room. To 
illustrate this, we now will determine the distribution of the busy period, the 
number of the customers, time of the first loss of the customer and the number 
of lost customers at time of the first loss. 

Corollary 5. Let P[6 = 1] = 1, bf.{x) = E [e~'^'^''^^^;br{x) < oo] be the Laplace 
transform of the busy period of type {r,x) of M^|G|1|B system. Then the fol- 
lowing relation is valid 

f,{s) + {l-f{s))S-l,_^.{x) 



b'{x) = '-^ ^ ^'^'^ ^ x>0, rGl,B + l, 

f{s) + {l-f{s))S% ' - ' 
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where S^{x) = X^^Lq^K^)' ^H^) ~ f^''" k < 0. The random variable br{x) 
is proper and 

Bbrix) = Er^a; - Brj + Brj [Sb - Ss-rix)] < oo, 
where Sk{x) = S^{x). 

These formulae were derived in [23]. To prove the corohary, it suffices to set 
A = in the equalities of Corollary [3j 

Corollary 6. Let P[5 = 1] = 1, lr{x) be the time of the first loss of the batch 
of the customers in the system M^|G|1|B, and lf.{x) = E [e~^'"'^^^ -Jrix) < oo] 
be the Laplace transform of lr{x). Then the following formula holds 

S ~\~ fJi \^ S ~\~ fJj J 

where r e 0,5 + 1, /c = S + l-r, = XliLV = ^IQfi+i-i- The random 

variable lr{x) is proper and has finite mathematical expectation. 

The function lf.{x) was found in [23] in a different form. 

Corollary 7. The distribution of the number of the customers in the system 
M'^IGIllB at time Ug is such that for r G 0,5 + 1, u £ 1,B + 1 

qUu)=A--^{s)+b^Ax) ^ , qUB + 1) = 1 '^'^''^ 



s + fj, — fib{s) ' s + n — fib{s) 



s + fj, — ^ b[s) 



Corollary 8. Let vrj = lim P[(i(.)(t) = i], i G 0,5 + 1 be the stationary 



t^oo 

distribution of the number of the customers in the system M^|G|1|B. Then 

-1 



1 + fiBr] ^ CLiQB-i , 

K i=0 / 



TTfi+l = 1 - VToQb, TTj = TToiQi - Qi-l), i € l,B, 

where the resolvent sequence of the process {Qk{x)}k(^z+ ^ Qk '= Qfc(O), Qo = 1 
is given by for s = 0. 

To prove Corollaries [HHSl it suffices to set A = in the equalities of Corollaries 
[2l[3]and Theorem [6l Note, that the formulae of Corollary [8] were obtained in [23]. 

Suppose that the system starts functioning from the state (r, x), r G 0, 5 + 1. 
Denote by ir,x the number of the lost customers at time of the first loss lr{x). 
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Corollary 9. The generating function Lf.^{z) = E [e ^^r(x)^ir,xj of the joint 
distribution {lrix),ir,x} is such that 



+ 



i=0 



Q%+i s + fi - nQ{s)/Q%^-^ 



k = B + l-r, (51) 



E 



— Slr{x) ■ — 



s 



^Axis) + ^{an+i - a„+i-i)^^ \s) 



1=1 



+ 



B+1 



Qb+1 S + fl- flQ{s)/Q%^^ 

where E'{z) = E [z^^*; x > , i € Z+. 



n G N, 



Proof. In accordance with the total probabihty law we can write the following 
system of equations for Lf.^{z) 



K,xiz) = V''{x,z,s) +Vr-^i{x,s)L^o^Q{z), r G 1,5 + 1 

B+l 



where 

V''{x,z,s) = E 



S + 



r=l 



(52) 



sx-.A^),T.^k =f>^(^^,,,s)-9MfB+i^o,z,s). (53) 



The equality <^ implies for the function f^{x,z,s) = E e'^'^'^^^'^z-^'^^^^; *B^(x) 
that 

k 



f'{x,z,s) = ^Y.^'('^ [At\s)-A'^,-^-Hs] 

1=0 

l-c(s) k{z)-k(c{s) 



+ QUx)F{c{s),z), (54) 
where F{c{s),z) = I'i^'^g^j^ '^^s-lt(c(s)' ■ Solving the system (j52]) and taking into 



account (j53j) for all r € 0, + 1, x > 0, we get 



B+l 



LlAz) = f^{x,z,s) + ii 



E^+\z)+ E aj^+i-(0,z,.)-p+i(0,z,s) 

Qkix) i=i 

Qb+i 



s + n- fiQ{s)/Q%^^ 



The first formula of the corollary follows from the latter equality and from 
Comparing the coefficients of z", n G N in both sides of ()5ip . we derive the 
second formula of the corollary. □ 
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